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ABSTRACT 

We present time-series spectroscopy of the rapidly oscillating Ap star HR 3831. This 
star has a dominant pulsation period of 11.7 minutes and a rotation period of 2.85 d. 
We have analysed 1400 intermediate-resolution spectra of the wavelength region 6100- 
7100A obtained over one week, using techniques similar to those we applied to another 
roAp star, a Cir. 

We confirm that the Ha velocity amplitude of HR 3831 is modulated with rotation 
phase. Such a modulation was predicted by the oblique pulsator model, and rules out 
the spotted pulsator model. However, further analysis of Ha and other lines reveal 
rotational modulations that cannot easily be explained using the oblique pulsator 
model. In particular, the phase of the pulsation as measured by the width of the Ha 
line varies with height in the line. 

The variation of the Ha bisector shows a very similar pattern to that observed in 
a Cir, which we have previously attributed to a radial node in the stellar atmosphere. 
However, the striking similarities between the two stars despite the much shorter 
period of a Cir (6.8 min) argues against this interpretation unless the structure of the 
atmosphere is somewhat different between the two stars. Alternatively, the bisector 
variation is a signature of the degree £ of the mode and not the overtone value n. 

High-resolution studies of the metal lines in roAp stars are needed to understand 
fully the form of the pulsation in the atmosphere. 
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1 INTRODUCTION 

The rapidly oscillating Ap (roAp) stars pulsate in high- 
overtone non-radial p-modes (Kurtz 1990; Martinez & Kurtz 
1995; Matthews 1991, 1997; Shibahashi 1991). The oscilla- 
tion spectra of many of these stars reveal frequency multi- 
plets - usually triplets - with equal frequency separations 
of the components. The separations are equal to, or very 
nearly equal to, the rotation frequency (determined from 
the variation of spectrum and light as the star rotates). 

One explanation of the observed triplets is that a set of 
rotationally perturbed m- modes is excited, e.g., i — 1 dipole 
modes with m = —1, 0, +1. The observed frequencies can 
be written as 



= v o - rn(l - C nt )v vt 



(1) 



where C n i is a constant that depends on the structure of 
the star (Ledoux 1951). For A-star models with the expected 
pulsation modes for roAp stars, Takata & Shibahashi (1995) 
found C n p ~ 0.003-0.010. However, in the best-observed 



roAp star HR 3831, Kurtz et al. (1992) were able to show 
that C n i < 2 x 10 -5 . This is two orders of magnitude less 
than the lowest theoretically predicted values, and suggests 
that the frequency splitting is precisely the rotational fre- 
quency. 

Another explanation for a frequency multiplet is a sin- 
gle pulsation mode (e.g., I = 1, m = 0) whose amplitude 
is modulated by rotation. This naturally produces a fre- 
quency splitting that is precisely equal to the rotation fre- 
quency. In some roAp stars, the amplitude of the pulsa- 
tion has been observed to modulate with rotation in phase 
with the magnetic variation. To explain this phenomena, 
the oblique pulsator model was proposed by Kurtz (1982), 
and has been extensively developed by Shibahashi & Takata 
(1993; Takata & Shibahashi 1994, 1995). In this model, the 
magnetic and pulsation axes are aligned but are oblique to 



In this paper, we use modulation to describe changes in am- 
plitude and phase with rotation, and variation for other types of 
changes. 
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the rotation axis. The observed amplitude modulation is 
then due to the variation in the angle between the pulsa- 
tion axis and the line of sight as the star rotates. 

In an alternative model, the spotted pulsator model 
(Mathys 1985), the pulsation and rotation axes are aligned 
(or there are radial I = modes), but the ratio of flux to 
radius variations varies over the surface because of differ- 
ences in the flux and temperature caused by spots associ- 
ated with the magnetic field. The observed radial-velocity 
(RV) amplitude is then constant, but the photometric am- 
plitude modulates as the star rotates. Therefore, the two 
models can be distinguished by RV observations. Matthews 
ct al. (1988) found RV amplitude variation in HR 1217 which 
favoured the oblique pulsator model, but this star is multi- 
periodic, so the observed variation could have been caused 
by beating among pulsation frequencies. Kurtz et al. (1994a) 
proposed that the two models could easily be distinguished 
by RV measurements of HR 3831. This roAp star has one 
dominant pulsation mode with a period of 11.7 min, has a 
large photometric amplitude modulation including a phase 
reversal (Kurtz et al. 1994a, 1997) and is bright enough for 
accurate RV measurements (V = 6.25). 

We obtained time-series spectroscopy of HR 3831 over 
one week in 1997 March. Preliminary analysis of the Ha RV 
measurements was presented by Baldry, Kurtz & Bedding 
(1998b; hereafter BKB) and showed that the radial- velocity 
amplitude is modulated with the rotation of the star. A fre- 
quency analysis showed a frequency triplet with the same 
spacing and amplitude ratios as contemporaneous photo- 
metric observations. This is expected in the oblique pulsator 
model, and in clear disagreement with the prediction of the 
spotted pulsator model. 

In this paper we describe more detailed analysis of our 
HR 3831 spectra using some of the techniques that were 
used by Baldry et al. (1998a, 1999; hereafter B98, B99) to 
study a Cir. These include cross-correlations of different 
wavelength regions, line-shifts of the Ha bisector and in- 
tensity measurements across the Ha line. While the oblique 
pulsator explains most of the observed rotational modula- 
tions, aspects of the spotted pulsator model and perturbed 
m-modes may also be required to explain our results. 



above 0.3mmag (Johnson B) and that are detectable in 
our data. These include a triplet around the principal mode 
= 1423.95 ^Hz, v = 1428.01 ^Hz, v +1 = 1432.07 ^Hz) 
and the first harmonic {2vq — 2856.02 /iHz). 

The splitting of the triplet (f ro t = 4.06 /iHz) gives the 
rotation period P ro t = 2.851976 d, assuming the oblique pul- 
sator model is correct. Unlike the triplet in a Cir (Kurtz et 
al. 1994b), where the amplitudes of the two side-lobes are 
about 10 percent of the principal amplitude, the amplitudes 
of f-i and are larger than the amplitude of the cen- 
tral frequency. This means that during the rotation cycle of 
HR 3831, the amplitude measured at uq goes through two 
maxima with two phase reversals. In terms of the oblique 
pulsator model, the star pulsates in an I — 1 dipole mode 
which is sufficiently oblique to the rotation axis that we see 
first one pole and then the other as the star rotates. 

Kurtz, Shibahashi & Goode (1990; Kurtz 1990) de- 
scribed a generalised oblique pulsator model in which the 
effects of both the magnetic field and rotation were taken 
into account (see also Dziembowski & Goode 1985, 1986). 
In this model, the perturbation to the star's eigenfrequencies 
by the magnetic field dominates, leading to the conclusion 
that the pulsation axis is locked to the magnetic axis. Two 
parameters related to the amplitude ratios in the frequency 
triplet were defined: 

P 1 = A + 1 + A - 1 =tanitan/3, (2) 
Ao 

where i is the inclination of the rotation axis to the line of 
sight, and (3 is the angle between the rotation axis and the 
pulsation axis; and 

p _ A+\ — A-l _ Cnt Vrot ,„s 

2 _ A +1 + A-! ~ v mmag _ v 0-)mag ' ^ > 

where the perturbation to the eigenfrequencies (by the mag- 
netic field) depends on \m\ such that 

From the measured photometric amplitudes for HR 3831 
covering 1993 to 1996 (Kurtz et al. 1997), we derive Pi = 
8.6 ± 0.2 and P 2 = -0.097 ± 0.004. 



1.1 Basic data for HR 3831 

HR 3831 (IM Vel, HD 83368, HIP 47145) has a binary com- 
panion with V — 9.09 with a separation of 3.29 arcsec. The 
visual magnitude of HR 3831 is often quoted as V = 6.17, 
but this is from the combined flux measurement (see Kurtz 
et al. 1994a). Using the correct value of V = 6.25, the Hip- 
parcos distance of 72.5 ± 4pc (ESA 1997), an effective tem- 
perature of 8000 ± 200 K (Kurtz et al. 1994a) and a bolo- 
metric correction of —0.13 (Schmidt-Kaler 1982), we obtain 

L = 13.4 ± 1.5L and R = 1.9 ± 0.2R Q .t 

The oscillation spectrum of HR 3831 has a well-known 
frequency septuplet around 1428 /iHz, plus frequency multi- 
plets at the harmonics (Kurtz et al. 1997). However, there 
are just four frequencies which have photometric amplitudes 

t In some papers (including BKB), the radius of HR 3831 is in- 
correctly quoted as 2.9R0. This may be due to propagation of a 
typographical error. 



2 OBSERVATIONS 
2.1 Spectroscopy 

We obtained 1400 intermediate-resolution spectra of 
HR 3831 using the coude spectrograph on the 74- inch Tele- 
scope at Mt. Stromlo Observatory (MSO), Australia. A log 
of the observations is shown in Table 1. The wavelength re- 
gion from 6100A to 7100A was observed with a resolution 
of about 1.5A and a dispersion of 0.49 A/pixel. The average 
number of photons/A in each spectrum was 270 000, with 
a minimum of 60 000. The exposure time was 100 seconds, 
with an over-head between exposures of 20 seconds. 

All CCD images were reduced to spectra using the same 
procedures that were used on the Mt. Stromlo data for B98. 
Several observables were defined from the spectrum (veloc- 
ities, intensities, bisector measurements), to produce a time 
series of measurements for each observable. For use in var- 
ious measurements, a template spectrum was defined from 
25 high-quality spectra (Fig. 1). 
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Table 1. Log of the spectroscopic observations of HR 3831 



UT-date 


No. of 


No. of 


Julian dates 




hours 


spectra 


-2450000 


1997 Mar 10 


1.37 


26 


517.88 - 


517.94 


1997 Mar 10 


3.61 


38 


518.12 - 


518.27 


1997 Mar 11 


10.03 


203 


518.88 - 


519.29 


1997 Mar 12 


9.59 


285 


519.88 - 


520.28 


1997 Mar 13 


2.22 


59 


520.88 - 


520.97 


1997 Mar 13 


1.10 


33 


521.27 - 


521.32 


1997 Mar 14 


10.57 


319 


521.87 - 


522.31 


1997 Mar 15 


10.40 


262 


522.88 - 


523.31 


1997 Mar 16 


9.40 


158 


523.87 - 


524.27 


1997 Mar 17 


0.53 


17 


524.87 - 


524.90 




6200 6400 6600 6800 7000 

Wavelength (A) 



Figure 1. The template spectrum of HR 3831 



2.2 Photometry 

For the last five years, astronomers at the South African 
Astronomical Observatory (SAAO) and the University of 
Cape Town (UCT) have been monitoring the oscillations in 
HR 3831 using high-speed photometry (Kurtz et al. 1997). 
For BKB, a subset of their data was selected which is centred 
on the time of the spectroscopic observations. These data 
were obtained on 26 nights (one hour of high-speed photom- 
etry per night) spanning the dates JD 2450402 to 2450618. 

The photometry was analysed to find the amplitudes 
and phases of the frequency triplet (see Table 2) in order to 
compare with the spectroscopic results. The amplitudes of 
the triplet are within 2a of the amplitudes from Kurtz et al. 
(1997). For further details, in particular showing the pulsa- 
tion amplitude and phase as a function of rotation phase, 
see BKB. 



Table 2. Photometric (Johnson B) amplitudes and phases" of 
the frequency triplet (BKB). From this data, Pi = 6.9 ± 0.7 and 
P 2 = -0.07 ±0.02. 



frequency 


amplitude 


phase" 


(jiHz) 


(mmag) 


(radians) 


1423.95 


1.91 ±0.05 


1.00 ±0.03 


1428.01 


0.52 ±0.05 


-1.16 ±0.10 


1432.07 


1.66 ±0.05 


1.03 ±0.03 



"The phases are shifted compared to the those given in BKB due 
to a different phase reference point (see Section 3). 
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Figure 2. R cw time series before and after high- pass filtering 



3 TIME-SERIES ANALYSIS 

To illustrate the various methods of analysing the time series 
for each spectroscopic observable, we look at an observable 
which has a high signal-to-noise ratio in the oscillation spec- 
trum. This observable, called R cw , is a ratio of the Ha core to 
wing intensity. In particular, we divided the intensity in a fil- 
ter with FWHM~ 4A by that in a filter with FWHM~ 31 A. 
This is similar to the R cw measurements used for B99 and 
BKB but using slightly narrower filters. 

Each time series from the spectroscopic measurements 
was high-pass filtered and a few outlying points were re- 
moved. The Rcw time series before and after high-pass fil- 
tering is shown in Figure 2. There is a variation in the raw 
data of about 1 percent during the night. This is possibly 
due to varying contamination of the HR 3831 spectrum from 
its binary companion. As the pair of stars rotates in the 
coude focal plane, there will be less contamination when the 
stars line up perpendicular to the slit and more contami- 
nation when the stars line up along the slit. Alternatively, 
the variation could be due to instrumental effects that are 
a function of elevation of the telescope. In either case, the 
change in the measured oscillation amplitude will be of the 
order of 1 percent over the night. This is not significant for 
our results, especially since there is no obvious systematic 
effect as a function of the rotation phase of HR 3831. 

Next, a simultaneous fit of the frequency triplet and 
the first harmonic was made to the high-pass filtered data, 
using a weighted least-squares fitting routine. The ampli- 
tudes and phases were fitted, using the well-known frequen- 
cies (see Section 1.1), by the function Asm{2irv(t — to) + 4>) 
where to = JD 2450522.51746. Our phase reference point to is 
equal to the reference point of Kurtz et al. (1997) plus 775 
times Prot- This maintains the same relationship between 
the phases of the frequency triplet. The fit of the frequency 
triplet to the R cw data is shown in Table 3 (the results for 
three sets of filters are shown). 

The R cw oscillation amplitude spectrum is shown in 
Figures 3-4. Figure 4 shows that the frequency triplet is re- 
solved, but with some aliasing or power shifting between the 
frequencies, as can be seen from the spectral window. The 
amplitude of the spectral window at ±^ ro t and at ±2^ rot is 
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Table 3. R cw amplitudes and phases of the frequency triplet in HR 3831. The observable 
R cw is the ratio between the mean intensity in a narrow filter with FWHM Fl and in a 
filter with FWHM F2 centred on Ho. For the three observables in this table, the average 
values arc about 0.65 (see Figure 2 for the raw data in the second case). 



Fl 

(A} 


F2 

(A} 


frequency 


amplitude 
(ppm) 


phase 
(radians) 


Pi 


P 2 


2.9 


23.5 


1423.95 
1428.01 
1432.07 


1340 ± 46 
172 ± 45 
1386 ± 46 


-0.82 ±0.03 
-1.28 ± 0.27 
-0.78 ± 0.03 


15.8 ±4.7 


+0.02 ±0.02 


3.9 


31.4 


1423.95 
1428.01 
1432.07 


1334 ± 41 
177 ± 40 
1139 ±41 


-0.71 ±0.03 
-1.58 ± 0.23 
-0.78 ± 0.04 


14.0 ±3.5 


-0.08 ±0.02 


5.9 a 


45. l a 


1423.95 
1428.01 
1432.07 


1107 ±41 
163 ± 41 
913 ±41 


-0.64 ±0.04 
-1.70 ± 0.25 
-0.78 ± 0.05 


12.4 ±3.4 


-0.10 ±0.03 



"This set of R cw results uses the same filters as the R cw results published in BKB. 
The results are slightly different due to a different weighting in the least-squares fitting 
routine, and there is also a change in phase reference point. 



about 20 percent. Also, the simultaneous fit gives amplitudes 
in ppm of 1330, 180 and 1140 for the triplet, whereas the 
amplitudes in the oscillation spectrum are 1100, 560 and 890 
(note that ratio between A_i and A +1 is about the same). 

Another way to consider the aliasing is that the data 
have incomplete and biased sampling of the rotation phase 
of HR 3831. A good diagnostic of the data is to plot the 
amplitude and phase of the pulsation (at i/o) as a function 
of rotation phase. This is better for interpretation assuming 
that there is one mode, that is amplitude modulated, rather 
than three distinct modes. We divided the data into 20 sep- 
arate time intervals between 0.5 and 3.6 hours long and in- 
cluding between 17 and 107 spectra. For each time interval, 
the amplitude and phase were measured at 1428.01 fiHz. The 
Rcw results are shown in Figure 5. The two lines represent 
different fits to the data: 

(i) The dashed line represents a modulation using the am- 
plitudes and phases of the R cw frequency triplet (from Ta- 
ble 3), which are from a six parameter fit to the complete 
time series. In other words, it shows the beating effect of 
the three close frequencies fo, f+i)- The good agree- 
ment between the modulation in amplitude and phase from 
the separate time intervals and the modulation from the fre- 
quency triplet supports the accuracy of the simultaneous fit 
to the triplet. 

(ii) The dotted line represents a fit which is obtained by 
scaling and phase shifting from the photometric analysis of 
the triplet (from Table 2). In effect, this is a two-parameter 
fit to the amplitudes and phases of the separate time inter- 
vals. The shape of the amplitude and phase modulation is 
the same as the beating effect of the photometric triplet, 
while one parameter is the scaling of the amplitude and the 
other parameter is the shift of the pulsation phase. In this 
case, a reasonable fit is obtained to the R cw data. 

In conclusion, we have clearly detected the triplet in 
the equivalent- width of Ha (using R C m) with the amplitudes 
of the three components in about the same ratio as seen 
in photometry. The R cw measurement, as in a Cir (B99), 



provides a high signal-to-noise spectral measurement of the 
pulsation. 



4 VELOCITIES OF WAVELENGTH BANDS 

In this section, we look at the velocity amplitude and phase 
of different wavelength bands using a cross-correlation tech- 
nique. The reduction method is the same as that used on 
a Cir, with a telluric band used as a velocity fiducial (Sec- 
tions 3.2-3.3 of B98). The spectrum of HR 3831 was divided 
into 90 bands, most having the same wavelength range as 
shown in Table 3 of B98, the difference being that the spec- 
trum of HR 3831 was taken from 6100A to 7100A, which is 
100A higher than the range used in the a Cir analysis. The 
results for 10 selected bands are shown in Table 4. 



4.1 Ha velocity 

First, we look again at the Ha velocity (previously analysed 
for BKB), as measured using the cross-correlation of band 
no. 87 (6545-6578A). This band has good signal-to-noise 
and should not be significantly affected by blending. The 
Ha velocity amplitude and phase as a function of rotation 
phase are shown in Figure 6. There is excellent agreement 
between the fit based on the measurement of the frequency 
triplet (Table 4) and the fit scaled and phase shifted from 
the photometric data (Table 2). This rules out the spotted 
pulsator model. In this model, the pulsation velocity ampli- 
tude and phase should be constant during the rotation of 
the star, which is clearly not the case. 

The excellent agreement between the Ha velocity and 
photometry, in terms of the relative modulation as a function 
of rotation phase, is reflected in the fact that the parame- 
ters Pi and P2 are nearly the same between the two time 
series (compare Table 2 with band 87 from Table 4) . Within 
the oblique pulsator model, these parameters are expected 
to be the same for different observables (see Section 1.1). If 
the frequency triplet were caused by three different modes, 
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Table 4. Velocity amplitude and phases for selected wavelength bands using a cross-correlation 
technique. The bands selected were those having a combined signal-to-noise a ratio of greater than 
nine for the frequency triplet, except only one band across the Ho line is included (band no. 87). A 
telluric band from 6864 to 688lA is used as a velocity reference (band no. 80). 



band 6 
range (A) 


frequency 
(Litiz) 


amplitude 
(ms ) 


phase 
(radians) 


Pi c 


P2 


Figure 


no. 13 
fii a n ft fti c^n 


1423.95 
1428.01 

1 /I *3Q fl7 
140Z.U I 


377 ± 51 
92 ± 50 

ZOU it 01 


0.82 ±0.13 
-0.67 ±0.57 

n so -U n on 
U.oy it U.zu 





-0.18 ±0.12 


— 


no. 14 

ft1 £0 O ft1 ft/1 ft 
DlOZ.y— 0104.0 


1423.95 
1428.01 

1 /I *3Q C\1 
140Z.U I 


524 ± 61 
26 ± 61 

07O _|_ CI 

Z i it 01 


2.45 ±0.12 

Ai _i_ n oq 

Z.41 it U.Zo 





-0.32 ±0.12 





no. 18 

fti Q/i n fti 07 ^ 
oiy4.u— oiy t .0 


1423.95 
1428.01 

140Z.U I 


1108 ± 69 

228 ± 68 

QOO -U ftO 

oyz it oy 


0.07 ±0.06 
-1.33 ±0.30 

n no xn is 
u.uy it U. lo 


6.6 ±2.3 


-0.48 ±0.07 


9 


no. 33 

OoZO.O— OoOO.Z 


1423.95 
1428.01 

140Z.U ( 


1140 ± 62 
222 ± 62 
ooo it oz 


0.54 ± 0.05 
-0.96 ±0.28 

n a ft J_ n 11 
U.4o it U. 11 


7.6 ±2.4 


-0.34 ±0.05 





no. 42 
ft/i 1 a n ft/i 00 

0414.U— 04ZZ.y 


1423.95 
1428.01 

140Z.U i 


600 ± 64 
118 ±64 

A -i- ft/1 
40o it 04 


1.26 ib 0.11 

-0.77 ib 0.57 
n 7o _i_ n 1 /i 

U. 1 U it U. 14 





-0.13 ±0.09 





no. 54 

OOZl.o— OOZo.Z 


1423.95 
1428.01 

140Z.U f 


348 ± 60 
88 ± 60 
a o*} -U ftn 

4Zv5 it OU 


1.64 ±0.17 
-0.44 ±0.75 
1 fto _i_ n 1 /i 

l.OZ it U. 14 





+0.10 ±0.11 





no. 58 

ftc:Qft Q fiftf|7 1 


1423.95 
1428.01 

±4:OZr.U < 


1378 ± 132 
131 ib 131 

Q/lft _|_ 1 qo 


2.36 ib 0.10 
-0.47 ib 1.49 

O QQ _|_ n 
z.oo ^ u.oy 


— 


-0.60 ±0.13 


— 


no. 81 d 
6881.5-6901.6 


1423.95 
1428.01 
1432.07 


63 ± 12 
16 ± 12 
68 ± 12 


-2.59 ±0.19 
1.99 ±0.80 
2.54 ±0.18 




±0.03 ±0.13 




no. 87 e 
6545.4-6578.2 


1423.95 
1428.01 
1432.07 


343 ± 22 
94 ±22 
309 ± 22 


1.58 ±0.06 
-0.81 ±0.24 
1.42 ±0.07 


6.9 ±1.8 


-0.05 ±0.05 


6 


no. 90 / 
7070.2-7079.5 


1423.95 
1428.01 
1432.07 


743 ± 87 
116 ± 87 
357 ± 87 


0.95 ±0.12 
-1.15 ±0.84 
1.28 ±0.25 




-0.35 ±0.12 





"The rms-noise is about 1.38 times the formal lcr uncertainty on the amplitude measurements. 
6 The band number relates to the bands used in B98. 

c Only the bands where the central amplitude is higher than 2<r are included. For the other bands, 
the uncertainty in this parameter is large and significantly non-Gaussian. 

d This band contains a significant number of telluric lines, see discussion in Section 5.1 of B98. 
e This band represents the velocity of the H« line. These results are different to those published in 
BKB because the band is wider and the telluric band used as a reference is different. 
'This is a new band not included in B98. 



then the ratios between the amplitudes of the modes (quan- 
tified by Pi and P2) would be expected to vary depending 
on the observable. Therefore, this agreement between the 
Ha velocity and photometry argues in favour of the oblique 
pulsator model rather than for different modes. 

However, this is a weak argument because the ampli- 
tude ratios could be nearly the same even in the case of dif- 
ferent modes. The main arguments against different modes 
come from the years of photometric analysis of HR 3831 
(e.g.: Kurtz et al. 1990, 1993, 1994a, 1997): (i) the frequency 
splitting is exactly or nearly exactly the rotation frequency 
(see Section 1), and (ii) the ratios between the amplitudes 
of the frequency triplet has remained nearly constant over 



time. For these reasons, it is considered unlikely that the 
observed frequency triplet in HR 3831 is caused by different 
modes, in particular rotationally perturbed m-modes. 

In BKB, a possible rotational phase lag was noted, be- 
tween the radial- velocity (RV) amplitude maximum and the 
photometric amplitude maximum, of about 0.06 rotation cy- 
cles. However, from the results shown in Figure 6, the RV 
maximum and photometric maximum are within 0.02 rota- 
tion cycles of each other. The change from the earlier results 
arises from the use of a different telluric reference wavelength 
region. In BKB, a large region from 6865 to 693lA was used, 
which included bands 80, 81 and 82, in order to maximise 
the signal-to-noise ratio. In this paper, we find a significant 
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Figure 3. R C w amplitude spectrum. The lower panel is the pre- 
whitcned amplitude spectrum after subtracting a simultaneous fit 
of the frequency triplet and the harmonic. 
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Figure 4. Close-up of Figure 3, with the dotted lines showing the 
frequencies of the triplet. The inset shows the amplitude spectrum 
of the window function on the same frequency scale, with the 
dotted lines having the same spacing as in the triplet. 



signal arising from the wavelength region of band 81 (see 
Table 4) and this is likely to be the cause of the phase lag 
noted in BKB. Band 80 alone is used as the telluric refer- 
ence band in this paper because it has the lowest percentage 
(about 2%) of absorption that is attributable to lines from 
the star (see discussion in Section 5.1 of B98). To summarise, 
the Ha RV amplitude is better represented by the results in 
this paper, whereas the results in BKB are slightly contam- 
inated by metal lines around band 81. 



4.2 Metal lines 

The velocity amplitude and phase in HR 3831 varied sig- 
nificantly between different wavelength bands, as was seen 
in a Cir (B98). Figures 7-8 show the amplitude versus pul- 
sation phase in HR 3831 for the 24 metal bands with the 
highest signal-to- noise ratios (plus one band across Ha). The 
results measured at the frequencies z/_i and v+i are plotted 
separately, while the amplitude and phase of vq are not plot- 
ted because the signal-to-noise ratio is significantly lower at 
this frequency. There are two aspects to the results which we 
will consider: the modulation of the measured pulsation with 
rotation and the variation of the amplitude and phase be- 
tween different wavelength bands. The first is characterised 
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Figure 5. R cw amplitude and phase of the central frequency 
as a function of rotation phase. The squares represent the data 
divided into 20 separate time intervals between 0.5 and 3.6 hours 
long. The vertical lines are error bars, while the horizontal lines 
show the length of the time intervals. The dashed line represents 
a fit based on the measurement of the frequency triplet from the 
complete time series, and the dotted line represents a fit which is 
scaled and phase shifted from the photometric frequency triplet. 
See Section 3 for further details of the time-series analysis. Note 
that the data are plotted twice. 



by the difference in amplitude and phase between v~\ and 
is+i. 

4-2.1 Rotational modulation 

There is a significant difference between the results for v-\ 
(Fig. 7) and v+i (Fig. 8). Firstly, the amplitudes vary be- 
tween and 1500 ms -1 for f_i, whereas the amplitudes 
for v+i are all below 800 ms -1 . Secondly, the phases range 
across 2n for V-i, while the phases for v+i range between 
and 7T except for one low S/N band. This is not as expected 
from the oblique pulsator model. In this model, for each 
band, the phase of V-\ and v + \ should be the same and the 
amplitudes should be nearly the same (A + i w 0.87 x A-i 
from Table 2). Therefore, these results suggest that the two 
frequencies measured are separate modes, rather than rota- 
tional 'side lobes' of one principal mode. However, there is 
strong evidence for the oblique pulsator model as described 
in Section 4.1. Therefore, we need to consider another expla- 
nation within the framework of the oblique pulsator model. 
We first look at some selected bands, in terms of their mod- 
ulation of amplitude and phase with rotation. 

Recall that the oblique pulsator model involves param- 
eters Pi and P2 (Eqns. 2 & 3). All the measured spectro- 
scopic observables in this paper produce values for Pi that 
are within 3<r of the photometric value. However, there are 
large uncertainties in this parameter because it depends crit- 
ically on the amplitude measured at uq, which is typically 
less than two or three times the noise level in our data. 

The parameter P2 can be measured with higher accu- 
racy since it does not depend on the small central amplitude, 
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Figure 6. Ha velocity amplitude and phase of the central fre- 
quency as a function of rotation phase (band no. 87). Line styles 
have the same meanings as in Figure 5. There is excellent agree- 
ment between the two fits, which rules out the spotted pulsator 
model. 
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Figure 7. Velocity amplitudes and phases measured at the fre- 
quency in HR 3831 for different wavelength bands. Bands 
with S/N greater than 3 are plotted with solid lines, those with 
lower S/N have dotted lines. The same 25 bands are plotted in 
this figure and in Figure 8, chosen so that the combined S/N is 
greater than 5. Only one band across Ha is included and is shown 
using a diamond. Note that the data are plotted twice for clarity. 



and we find it does vary significantly between different ob- 
servables. From the measured velocity amplitudes, there are 
three spectral bands (18, 33 and 58) that have values of P2 
that are formally 4-5a different from the photometric value 
(see Table 4). To test whether this is a significant change 
or just a product of aliasing, we looked at the modulation 
in velocity amplitude and phase as a function of rotation 
phase. Figure 9 shows this modulation for band 18. While 
the amplitude modulation could be fitted reasonably well by 
scaling from the photometry (dotted line), the discrepancy 
is obvious with the phase modulation. There is a noticeable 
pulsation phase change from rotation phase 0.8 through 1.0 
to 0.2, while the photometric phase is nearly constant. Sim- 
ilar effects are evident from the other bands, 33 and 58, for 
which the P2 value is significantly different from the photo- 
metric value. 

Figure 10 shows a plot of P2 versus the average am- 
plitude of frequencies and v+i. There is a noticeable 
tendency for higher amplitude bands to have lower values of 
P2 . This could be related to the line-formation depth or po- 
sition on the surface (i.e., spots). Alternatively, a systematic 
effect (e.g., blending) could be increasing the measured am- 
plitude at V-\ only causing an increase in average amplitude 
and a decrease in P2. 

As well as the parameters Pi and P2, which relate the 
amplitudes of the frequency triplet, there are also phase dif- 
ferences between the frequencies. For many bands, there is 
significant difference between the phases measured at v-\ 
and v+i. In our measurements, the equality of phase 0_i 
and ci+i means that the amplitude maxima occur at ro- 
tation phases 0.0 and 0.5 with phase jumps in between. 
Therefore, any phase difference, <f>d — <fr+i — <p-i, implies 
that there is a shift of the amplitude maxima. For example, 
there is a small shift in the amplitude maxima of band 42 
(<f>d = -0.56 ± 0.18) and larger shifts, but with lower S/N 
in the amplitude spectra, of bands 48 (4>d = —2.29 ± 0.31) 
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Figure 8. Same as Figure 7, but for the frequency v+i. 



and 82 (4> d = 2.99 ± 0.25). Figure 11 shows the modula- 
tion for band 48. Note that the pulsation-phase jumps occur 
out of (rotation) phase with the photometric pulsation-phase 
jumps. 

These results are not sufficient evidence to discard the 
oblique pulsator model, because they could be explained by 
variations in the ratios between the metal lines in the bands, 
due to spots, as the star rotates. Alternatively, there is a 
possibility that there is an actual change in the pulsation 
phase and amplitude (as opposed to a pseudo-change caused 
by a variation in blending) associated with spots on the 
star. So while the spotted pulsator model is ruled out in 
its original form, spots may be having an influence on the 
pulsation. 
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Figure 9. Velocity amplitude and phase of band no. 18 as a 
function of rotation phase. Line styles have the same meanings as 
in Figure 5. 
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Figure 11. Velocity amplitude and phase of band no. 48 as a 
function of rotation phase. Line styles have the same meanings as 
in Figure 5. 
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Figure 10. Plot of the parameter P2 versus the average ampli- 
tude of frequencies and with lcr-error bars, for the 25 
bands that are shown in Figures 7 and 8. Points with solid lines 
represent the higher S/N bands given in Table 4. Note that the 
four lowest amplitude bands in this figure are contaminated by 
telluric lines. For these bands, the cross-correlation measurement 
of the velocity amplitude is reduced compared to an uncontami- 
nated band (with equivalent stellar behaviour). 
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Figure 12. The Ha line in HR 3831. The dotted lines divide the 
25 contiguous sections used in the bisector-velocity and width 
analysis (Section 5). 



nation that includes running waves is necessary. Alternative 
explanations could involve spots on the surface or blending, 
which causes the measured phase to deviate from the true 
velocity phase (for further discussion see Section 6.1.2). 



4-2.2 Amplitude and phase variations between bands 

In B98, it was suggested that the amplitude and phase varia- 
tions between wavelength bands in a Cir could be explained 
by a radial node in the atmosphere. We see amplitude and 
phase variations in HR 3831 that are somewhat similar in 
the case of (Fig. 7), and with smaller but still significant 
variations in the case of v+\ (Fig. 8). However, there is no 
obvious division between two sets of bands that are pulsating 
in anti-phase with each other, as was seen in a Cir (B98) . If 
these HR 3831 results represent true velocity amplitudes and 
phases at various depths in the atmosphere, then an expla- 



5 Ha PROFILE VARIATIONS 

In a Cir we found intriguing behaviour in the Ha line that 
implied the existence of a radial node in the atmosphere 
(B99). We now describe a similar analysis for HR 3831. The 
Ha line in each spectrum was divided into 25 contiguous 
horizontal sections, between relative intensities 0.35 and 0.9 
(Fig. 12). For each section, a bisector velocity (average po- 
sition of the two sides) and a width were calculated. For the 
velocity measurements, a telluric band (no. 80) was used as 
a fiducial. 
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Figure 13. Amplitudes and phases of the pulsation for the bi- 
sector velocity at different heights in the Ha line. Squares with 
solid lines represent v—\ and diamonds with dotted lines repre- 
sent For each measurement, the vertical line is an error bar 
while the horizontal line shows the extent of the section in the 
Ha line. 

5.1 Bisector velocities 

The velocity amplitude and phase of the pulsation as a func- 
tion of relative intensity in the Ha line are shown in Fig- 
ure 13, with measurements at frequencies V-\ and v+\. The 
amplitudes and phases are in good agreement between the 
two frequencies, which means that the shape of the rota- 
tional modulation is similar for each velocity measurement. 
Figure 14 shows the rotational modulation from the veloc- 
ity at a height of 0.40. It is similar to the modulation from 
the cross-correlation velocity of Ha (Fig. 6), but with about 
twice the amplitude. This is expected because the cross- 
correlation velocity is a weighted average of the bisector ve- 
locities and the velocity amplitude at 0.40 is higher than 
average. 

The variation in amplitude and phase as a function of 
relative intensity (Fig. 13) is similar to that seen in a Cir, 
particularly the drop in amplitude between heights 0.4 and 
0.6 (a comparison is shown in Figure 15). However, assuming 
that the bisector line-shift represents the velocity at different 
heights in the atmosphere, we expect to see less relative 
variation in the bisector-velocity amplitude in HR 3831 than 
a Cir. This is because the separation between radial nodes 
in the atmosphere is larger in HR 3831 due to the fact that 
it pulsates at a lower frequency and that the two stars have 
similar fundamental properties. 

Why are the behaviours of a Cir and HR 3831 so simi- 
lar? Two possible explanations are: (i) the Ha line forming 
region in HR 3831 is more extended than in a Cir and co- 
incidently the Ha bisector behaves in a similar way; (ii) the 
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Figure 14. Velocity amplitude and phase of the Ha bisector at 
height 0.40 as a function of rotation phase. Line styles have the 
same meanings as in Figure 5. 
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Figure 15. Comparison of the Ha bisector variations between 
HR 3831 and a Cir. This figure is essentially a normalised over- 
plot of Figure 13 of this paper and Figure 7 of B98, except that for 
HR 3831 the measurements for V—\ and v+\ were averaged and for 
a Cir the Mt. Stromlo and La Silla data were combined. Addition- 
ally, the amplitudes were normalised by dividing by an average 
taken between heights 0.40 and 0.45 (519 ms -1 for HR 3831 and 
267 ms -1 for a Cir) and the phases were shifted. 
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Figure 16. Amplitudes and phases of the pulsation for the width 
at different heights in the Ha line. Points with solid lines represent 
v—\ and points with dotted lines represent the v+\. For each 
measurement, the vertical line is an error bar while the horizontal 
line shows the extent of the section in the Ha line. The dashed line 
shows the effect on the width amplitude of an oscillation with an 
EW amplitude of 1000 ppm, with the profile variation described 
in Section 4.2 of B99. Note that 10 mA is equivalent to a velocity 
of 457 ms^ 1 at Ha. 

assumption that the bisector represents velocities at differ- 
ent heights is incorrect. 

The variation in amplitude and phase is less similar to 
a Cir above a height of 0.7. For instance, there is a notice- 
able amplitude peak and phase reversal around 0.73 in the 
bisector of a Cir which is not evident in HR 3831. The bisec- 
tor measurements above a height of 0.7 in a Cir are partially 
affected by blending (Section 3.2 of B99). Therefore, some 
of the differences in the bisector- velocity results between the 
two stars are probably due to a difference in blending. 

To summarise, there is no significant difference in the 
rotational modulation of the bisector velocity between var- 
ious heights. For the measurements below a height of 0.7, 
the relative variation in amplitude and phase is very similar 
to that observed in a Cir, which is puzzling. At each height, 
the amplitudes at the two frequencies in HR 3831 are about 
twice the amplitude of the principal frequency in a Cir. 

5.2 Widths 

Changes in the width of the Ha line in HR 3831 at vari- 
ous heights are shown in Figure 16. The width amplitude 
and phase represent oscillations about a mean width, with 
measurements at two frequencies. These measurements are 
related to intensity and EW measurements, with the line re- 
solved vertically rather than horizontally. The advantage of 
width measurements is that the width is naturally indepen- 
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Figure 17. Width amplitude and phase of Ha at height 0.36 as 
a function of rotation phase. Line styles have the same meanings 
as in Figure 5. 

dent of the bisector velocity, i.e., the width is the difference 
in position between two sides of a line whereas the veloc- 
ity is the average position. In order to compare with EW 
measurements, the dashed line shows the simulated width 
amplitudes for the profile variation described in Section 4.2 
of B99 (see also Figure 6 of that paper) with an EW ampli- 
tude of 1000 ppm. 

The most striking feature of Figure 16 is the differ- 
ence in the measurements between the two frequencies. This 
means that, unlike for the bisector velocities, the shape 
of the rotational modulation is varying between different 
heights. We show three examples: 

(i) Figure 17 shows the modulation in width amplitude 
and phase at a height of 0.36. The pulsation width phase 
increases steadily during the rotation and P2 is large (0.54 ± 
0.05). 

(ii) Figure 18 shows the modulation at height 0.42. There 
is good agreement with the fit derived from the photometric 
triplet and P2 is —0.03 ± 0.04, which is close to the photo- 
metric value. 

(iii) Figure 19 shows the modulation at height 0.53. The 
width phase decreases steadily during the rotation and P2 
is -0.63 ± 0.07. 

To check the lowest measurement at height 0.36, which 
is sensitive to intensity changes in the core of the line, 
we took an intensity measurement using a filter with a 
full-width half-maximum of about lA. This measurement 
showed a similar rotational modulation, with the parameter 
P 2 = 0.33 ± 0.03. 

Something unexpected is certainly happening here. The 
Ha width changes at the two extremes, heights 0.36 and 
0.53, are in phase at rotation phases around 0.5 and 1.0, 
but in anti-phase at phases around 0.25 and 0.75. At both 
heights, the relative amplitude modulation (as defined by 
the dashed line) is less than for other observables. 

Looking at Figure 16, the simple explanation is that we 
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Figure 18. Width amplitude and phase of Ho at height 0.42 as 
a function of rotation phase. Line styles have the same meanings 
as in Figure 5. 
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Figure 20. Amplitudes and phases of the pulsation in 7 Equ for 
different metal lines, using results from Kanaan & Hatzes (1998). 
Lines with A/a a greater than 4 are plotted with solid lines, those 
with 2.5 < A/cr a < 4 are plotted with dotted lines. Note that the 
data are plotted twice for comparison with Figures 7 and 10 of 
B98 and Figures 7 and 8 of this paper. 
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Figure 19. Width amplitude and phase of Ha at height 0.53 as 
a function of rotation phase. Line styles have the same meanings 
as in Figure 5. 

are seeing two modes with different Ha profile variations. 
An alternative explanation - if HR 3831 is pulsating in an 
oblique mode - is that, as the star rotates, the profile varia- 
tion modulates due to viewing different aspects of the same 
mode. 

6 DISCUSSION AND CONCLUSIONS 
6.1 Velocity amplitudes and phases 

6.1.1 Another case: 7 Equ 

Kanaan & Hatzes (1998) observed radial-velocity variations 
© 2000 RAS, MNRAS 000, 1-13 



in the roAp star 7 Equ. Their observations consisted of high- 
resolution spectra in the range 5000-6000A, covering 3.5 h 
on 1994 September 13. There are four modes in the range 
1330-1430 /iHz (Martinez et al. 1996), which are unresolved 
in that data set. Kanaan & Hatzes fitted a sine curve of fixed 
frequency (1380 piHz), but variable amplitude and phase, to 
the velocity measurements of each individual line, effectively 
treating the modes as a single pulsation. They discovered 
that the velocity amplitude varied significantly from line to 
line (see B98 and Section 4 for similar results on a Cir and 
HR 3831). 

We analysed some of these results on 7 Equ, given in 
Table 1 of Kanaan & Hatzes (1998), in order to compare 
with the results on a Cir and HR 3831. From their table, 
which lists 70 metal lines, we took the amplitude A (m/s), 
a a (m/s) and the 'Time of Maximum' (JD). Note that the 
times of maximum quoted by Kanaan & Hatzes cover sev- 
eral pulsation cycles. We converted the times to phases us- 
ing the pulsation period of 0.008387 d, and the phases were 
shifted so that a weighted mean was approximately zero. The 
phase uncertainty was taken to be the arcsin (a a /A). This 
phase uncertainty is in good agreement with the a b (days) 
quoted in their table for lines where A/a a is greater than 
2. Figure 20 shows amplitude versus phase for the 19 lines 
with the highest signal-to-noise ratio. The higher amplitude 
bands, plotted with solid lines, cover a range of 0.8 radians in 
pulsation phase (—0.26 to +0.54 in the figure). This is signif- 
icantly less than the variation in phase seen in a Cir (Fig. 7 
of B98) and HR 3831 (Fig. 7), where the bands cover virtu- 
ally the whole range of phases (2ir). Note that the rms-noise 
level in an oscillation spectrum is about 1.38a a , therefore, 
A/a > 4 is comparable to A/noise > 3. 

The difference between the results for 7 Equ (Kanaan 
& Hatzes) and those for the other two roAp stars, in terms 
of the variation in phase between different metal bands, 
could be due to: (i) wavelength region studied, 5000-6000 
vs. 6000-7000, (ii) spectrograph resolution, 0.2A vs. 1.5A, 
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(iii) pulsational characteristics of the stars. Further high- 
resolution studies of these roAp stars will be able to distin- 
guish among these causes. 

6.1.2 Depth and surface effects 

What are the causes of the velocity amplitude and phase 
variations, and of the Ha bisector variations? 

The main hypothesis is that there is a significant change 
in pulsation amplitude and phase with geometric depth in 
the atmosphere of an roAp star, with a radial node of a 
standing wave situated in the observable atmosphere. Ab- 
sorption lines are formed at different depths and therefore 
sample different parts of the standing wave. This can ex- 
plain a complete range of observed amplitudes plus a phase 
reversal. 

An alternative hypothesis for these variations is that 
the amplitude and phase of a non-radial pulsation mode will 
vary over the surface of a star, and that each spectroscopic 
observable is related to an integral over the surface. There- 
fore, if the integrals are sufficiently different between observ- 
ables, the measured amplitude and phase may vary signifi- 
cantly. These variations could be related to limb- darkening 
and / or spots. 

Neither hypothesis can fully explain the range of phases 
seen in a Cir and HR 3831. Blending effects may cause de- 
viations from the true velocity phases that are in addition 
to any atmospheric depth and surface effects. The smaller 
phase variations seen in the high-resolution results of 7 Equ 
lend some support to this. Blending can be considered as 
blurring the distinction between velocity and temperature 
changes in a star. 

Surface effects would be a good explanation if roAp 
stars were pulsating in modes with £ > 3 because, for these 
modes, there are several patches on the surface that have 
alternating phase. However, the frequency triplets in many 
roAp stars suggest that dipole modes (i.e., I = 1) are pre- 
dominant. For these modes, the measured velocity ampli- 
tude and phase are not expected to vary significantly be- 
tween different observables. For example, for a dipole mode 
with an axis that is aligned near to our line of sight, the 
pulsation phase will be the same across most of the viewed 
surface. 

The main arguments for depth effects being the prin- 
cipal cause of the amplitude and phase variations are: (i) a 
radial node is plausibly situated in the observable atmo- 
sphere (Gautschy, Saio & Harzenmoser 1998), and (ii) the 
bisector velocity varies with height in the Ha line. In Sec- 
tion 3.3 of B99, the possibility of systematic errors causing 
the bisector-velocity variations was discussed. Here, we make 
a simple calculation to test the accuracy of these velocities 
below a relative intensity of 0.7 in the line. 

For a Cir and HR 3831 (f-i and v+i), we averaged 
the bisector-velocity amplitudes below a height of 0.7 (giv- 
ing half weight to the amplitudes between 0.6 and 0.7) in 
order to compare with the cross-correlation measurements. 
This is a plausible estimate for the cross-correlation because 
the slope of the Ha profile is steep and nearly constant be- 
tween 0.4 and 0.6, and is less steep above 0.6. For a Cir, we 
obtained an estimate of ~170ms~ 1 while the actual cross- 
correlation measurements gave 168-182 ms _1 (depending on 
the band, nos. 85-88). For HR 3831, we obtained estimates 



of ~340 and ~290ms _1 for the two frequencies while the 
actual measurements gave 338-361 and 296-318 ms -1 re- 
spectively. The good agreement, between the estimates ob- 
tained from the bisector analysis and the velocity amplitudes 
from the cross-correlation measurements (using four differ- 
ent bands), argues that the bisector velocities are accurate. 

However, there is a similarity between the a Cir and 
HR 3831 bisector results (see Figure 15) that is possibly in- 
consistent with depth effects being the cause of the bisector 
variations. This is because the pulsations of these two stars 
have significantly different periods (7 and 12 minutes) and, 
therefore, have different vertical wavelengths, assuming that 
the sound speed is approximately the same between their 
atmospheres. If depth effects are causing the bulk of the 
observed bisector variations, the Ha line is formed over a 
larger vertical distance, and/or the sound speed is lower, in 
the atmosphere of HR 3831. 

Again, surface effects could explain the bisector varia- 
tion in a Cir and HR 3831 if they were pulsating in modes 
with I > 3 (see Hatzes 1996). Therefore, we are left with 
two relatively straightforward explanations for the ampli- 
tude and phase variations: (i) the stars pulsate in dipole 
modes and the variations are caused by depth effects, or 
(ii) they pulsate in I — 3 modes and the variations are caused 
by surface effects. A more complicated explanation could in- 
volve a combination of standing waves, running waves, sur- 
face effects, distorted modes and blending. 

6.2 Rotational modulation in HR 3831 

There are many examples of rotational modulation shown in 
this paper. The large variation in the shape of the modula- 
tion between different observables was an unexpected result, 
for which we see three possible explanations: 

(i) The oblique pulsator model is wrong and the observed 
frequency triplet is caused by three different modes. 

(ii) The oblique pulsator model is correct and the varia- 
tions are caused by systematic errors in the measurements. 

(iii) The oblique pulsator model is correct but needs to 
be modified to include effects due to spots or some other 
mechanism. 

In the first case, the relative amplitude between dif- 
ferent modes would be expected to vary for different ob- 
servables. This is the natural explanation for the differences 
between V-\ and v + \ of the metal- line velocity amplitudes 
and phases (Figs. 7 and 8) and of the Ha-width amplitudes 
and phases (Fig. 16). However, it would require some ro- 
tational phase-locking process to account for the frequency 
separations being equal to the rotation frequency. 

In the second case, no simple systematic error could ex- 
plain the variation in the shape of the rotational modulation. 
Any such error would need to be varying with the rotation. 
For instance, a variation in blending could possibly account 
for the rotational modulation of the metal-line amplitudes 
and phases, but it could not account for the Ha-width vari- 
ations. 

For the third case, we consider the Ha line because it 
is the least affected by blending. First, note that the in- 
tensity modulation, using a filter with a FWHM of about 
6A, is fitted well by the rotational modulation derived from 
the photometric triplet. This means that the total intensity 
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variations below a height of about 0.62 are in agreement 
with the oblique pulsator model. The bisector velocities are 
also in agreement, so it is only the width variations below 
0.62 that are not in agreement with the model. Perhaps this 
could be explained by hot spots in the upper atmosphere 
that cause significant amplitude and phase differences be- 
tween pulsational temperature changes in and around the 
spots. As the star rotates, the effect on the Ha profile would 
vary, leading to a hybrid model combining spots with an 
oblique pulsator. 



6.3 Future work 

It is clear that the high-resolution study (R<;40000) of differ- 
ent metal lines will be necessary to improve the understand- 
ing of the pulsation and the structure of the atmosphere in 
roAp stars, because of the problems associated with blend- 
ing at lower resolution (R ~ 5000). ft will then be possible 
to compare the velocity amplitudes and phases between dif- 
ferent atoms and ionization states. Theoretical work needs 
to be done to calculate formation depths for individual lines, 
including the effects of diffusion and the magnetic field. The 
aim is to build up a coherent picture of the pulsation in the 
atmosphere. 

Of the results on a Cir and HR 3831, the Ha pro- 
file measurements are least affected by blending and cannot 
easily be improved using high-resolution spectroscopy. The 
profile measurements require a good continuum fit across 
the Ha line, which is about fOOA wide in these A-stars. 
For this reason, a spectrum which is stable across at least 
400A is preferable. This is harder to obtain with high- 
resolution spectroscopy due to the blaze pattern and the size 
of CCD detectors. Although improvement could be made 
on the results using intermediate-resolution spectroscopy 
(R ^15000), modelling the formation of the Ha line is more 
important. Determining depth and surface effects on the Ha 
profile will complement high-resolution spectroscopic and 
theoretical studies of metal lines. 

Photometric oscillation spectra could be considerably 
improved from observations of roAp stars using small space 
telescopes. Two such missions plan to include roAp stars 
as part of a project to measure oscillations in nearby stars 
(MON$, see Baldry 1998; and MOST§). If successful, it will 
result in the detection of new modes and improved astero- 
seismology of the brightest roAp stars. 

The spectroscopic study of the pulsations in roAp stars 
has produced more questions than answers, due to the com- 
plexity of the observed spectral changes. For the same rea- 
son, the 'amount of information' obtainable is very large 
(e.g., variations in amplitudes and phases) compared to 
other pulsating stars. Therefore, roAp stars are prime stel- 
lar objects for testing diffusion, magnetic field and pulsation 
theories. 
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